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Transport coG“^r ic i noneltj therncl conductivity^ viccosity 

and diffusion constant, of external field-free fully ionized plasr-.as 
ere exui. lined in, the f rat.ie'rorli of a . inetic description oc zed on .'iJC/iT- 
Iiierarrliy o/ cc>:otionso 

Due to long 'ranc ,3 nature of the Couloizb jorce, leading to 
nultiplo collisions and interpe rticlc-corrclc t ions , a closed forr>. 
hinetic description oj plasiia : com only oe obtai '2d by nafiing a series 
0/ nimpli f yinn assumptions^ Sev^'^rol such schchies reported in literature 
are rev iev^edo 

It IS found that the correct functional dependence of placraa 
transport coefficients on plasma properties, vizo tepip^rature, dcncityp 
charge ana mass of the constituents, con be obiainrd by ixahing certain 
verij siiiiple assuhiptions renen ding the single , article and tizO’ par tide 
dijtribviion functions, loithin the fratccuo2''h of hydrozlynaaio ap.ro::i- 
niotionu Absoluhn nuiicrical values can then be calculated by co-tperizon 
unlh iUoasurchionta or hiore advanced Lheorvo 
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INTRODUCTIG?!^ DAChCJlQUhO AND OUTLINE 
1^1 Introduction 

A plasma is a quasineutral collection of positively 
and negatively charged particles. In addition to charged 
particles, neutral particles may also he present, in -which 
case the plasma is called partially ionized. If no neutral 
particles are present the plasma is called fully ionized. 

The transport processes in a fully ionized plasma are the 
subject of study in the present work. 

For a fully ionized plasma, the ionization-recombina- 
tion behaviour and collision with neutral atoms need not be 
modelled. In our work we will also ignore the relativistic 
and gravitational effects, and take plasma free from external 
electric and magnetic fields. 

Plasma behaviour can be studied by employing one of 
the following models ; 

1. Individual particle model, 

2. Kinetic model, 
and 3. Fluid model. 

Even with the availability of extremely fast computers, the 
individual particle model can not be used in almost all cases 
of practical interest because an astound%gly large number of 
equations are involved in order to study the system. Thus 
one has to rely on the kinetic or fluid models. 
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In the fluid or hydrodynamic model, macroscopic 
conservation equations for mass, momentum, energy etc are 
derived phenomenologically without any reference to the macro- 
scopic distribution function. The resulting equations, 
however, can not be solved without making additional assump- 
tions. The most Important of these assumptions is the intro- 
duction of ''linear rate equations" into the hydrodynamic 
equations. A potential gradient is associtited with the flux 
or net flow rate of a property based on experimental evidence 
and experience, and a linear relationship between the flux 
and the potential gradient is assumed. The transfer process 
takes place in the direction of decreasing potential. These 
phenomenological rate equations are s Pick's equation of 
diffusion, Fourier's equation of heat condi:iction, Fewton's 
equation of viscosity and Ohm's equation of electrical conduc- 
tion. 

This method, however, has one serious limitation. The 
coefficients in the rate equations, namely, the diffusion 
coefficient, thermal conductivity, the coefficient of viscosity 
and the electrical conductivity, are to be determined purely 
by experiment. That is, this approach does not provide any 
theoretical values of the transport coefficients, nor does it 
give the functional relationship between the transport coeffi™ 
cients and other macroscopic properties of the system. 
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In the Kinetic model, however, one derives the hydro- 
dynamic equations by taking appropriately weighted moments of 
the kinetic equation governing the microscopic behaviour of 
the system. When the kinetic equation is solved in the hydro- 
dynamic approximation, theoretical values ox xhe transport 
coefficients can be obtained. Also rate equations do not 
have to be introduced adhoc, as in the phenomenological method. 

The information obtained about the transport coeffi- 
cients from the kinetic theory, however, depends on the 
assumptions made in solving the kinetic equation. It is shown 
in this work that for small deviations from equilibriijm tne 
functional dependence of the transport coefficients of fully 
ionized plasmas can be deduced in a relatively simple manner. 
Absolute numerical values can then be calculated by comparison 
with measurements or more advanced theory. 

1 c 2 Background and Lite rature^Reulm^ 

The mathematical theory of thermal conduction, viscosity 
and diffusion in neutral gases till recently has been based 
almost entirely on Boltzmann kinetic equation [l]. The equation 
is applicable only to dilute neutral gases because of the 
assumptions made in deriving the equation, namely, 

(a) Interparticle forces are short-range. A consequence 
of this is that the mean free path between collisions is much 
larger than the range of interparticle forces, or the mean free 
time is much larger than the collision time. 
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(To) The collisions -.re only binary. 

(c) The particles are statistically indepenc* ”>t. 

Several methods for solving the Boltzmann equation 
have been developed, for example those of Chapman -Enskog, 
Burnett and Grad [2]. Enskog has also applied Boltzmann 
equation to dense gases. The Boltzmann equation rs, however, 
completely unsuitable for plasmas, because none of the above 
mentioned assumptions are valid due to the long range nature 
of the Coulomb force. In particular, the following effects 
have to be explicitj.y accoimted for in the case of plasmas; 

(a) The particles are no longer statistically 
independant . In other words they are statistically correlate^. 

(b) All collision in reality are multiple collisions. 
In fact, in some treatments oi plasmas, for example, Bokker- 
Planck theory, short range binary collisions can be ignored. 

(c) No well-defined mean free path or time exists for 
plasmas because the particles are interacting all the time. 

Bogoliubov [ 3 ] has presented a. ne';'r approach to the 
formulation of kinetic theory based on liouville equation [ 4 ] 
for the N~particle distribution function. By defining the 
reduced s-particle distribution functions, it is shown that 
the liouville equation can be put in the form of the completely 
equivalent system of equations, called BBGKT hierarchy of 
equations [5,6]. It will be discussed in detail in ''haptei 2. 
By making certain assumptions about the time evolution of the 
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distribution functions and about the nature of the correlations 
in the past, Bogoliubov has shovrn that, as a first approzimation, 
the Boltzmann equation can be derived from the ELGEJ equations. 
Choh and Uhlenbeck [2] have solved the BSGirf equations in hi^er 
order approximations, thus accounting for the multiple collisions 
and correlation effects for dense gases. 

Prieman and Sandri [8,2,10] have also developed a theory 
to solve the BBGKI equations with somewhat less restrictive 
assumptions, particula.rly, about the correlations in the pa,st, 
Rostoker, Rosenbluth and Ischikawa have developed another 
approach for truncating the BBGKf hierarchy of equations [11,12, 
13]- This is based on the concept of cluster expa^nsion well- 
knox^n in statistical physics [14,15]. '4u [2] has shown that 

the essential features of all these approonhes are equivalent, 
eventhough, the mathonatical s’;eps involved are quite different. 

Guernsey and Lenard [ 2,17] have extended Bogoliubov ’s 
theory for neutral gases to plasmas. The essential difference 
involved is in the expansion parameter used. For the case of 
noutral gases the density of gas is used as an ezpar ion para- 
meter, For plasmas, it is the ratio of potential energy to 
kinetic energy, cal3 ed tbe plasma parameter, which is used for 
expansion. This will be discussed in more defcail ir bapter 2. 
For small values of plasma parameter, these authors have worked 
out the equilibrium! one-particle and two-particle distribution 
fucnctions. These will be utilized 'in cur xfork in’'hapte:- 
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Sundaresan, Wu and Rosenburg [18,19? 20] liave applied 
Bogoliubov’s theory to inhomogeneous plasmas and derived an 
expression for thermal conductivity of fully ionized plasmas 
for small tenperature gradients* The dominojit term in their 
expression [20] agrees well with the expression derived earlier 
by Spitzer and Harm [21,22] on the basic of Fokker-Planck 
equation. 

In this work we have derived expressions for thermal 
conductivity, viscosity and diffusion coefficient of a fully 
ionized plasma by utilizing one-particle and two-particle 
equilibrium distribution functions as obtainod by other 
authors [2 ,17]. Small deviations about the equilibrium 
solutions are assuraed due to small gradients in temperature, 
density and macroscopic or hydrodynamic velocity. The wei^ted 
integrals of the kinetic equacion are obtained yielding 
appropriate expressions for the transport coefficients. The 
mathematics involved in this approach is considerably simpler 
than that needed in other approaches. 

lo3 Outline of t he Pre!jent_J/ork 

InC hapter 2, a detailed description of the kinetic 
theories of gases and particularly plasmas is giver The 
emphasis is on BB&ET approach. The equals brium form of the 
one-particle and two-particle distribution functions is 
presented and discussed, InC hapter 3, the hydrodynamic 
equations are derived by taking the weighted moments of the 



kinetic equation, and the transport coefficients are identi- 
fied. In Chapter 4, the expressions for the transport coeffi- 
cients are derived based on simplifying assumptions mentioned 
above. Chapter 5 presents a summary of this work, reiterating 
the main conclusions and suggesting areas that need further 
work. 
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2IU!-, TIC OJhCnif TIOiJ 


'^JTPAL 


.'IS 


'D PL/.Ij, /15 


° " -^SSiypll-Boltzimnn Idnetic Theory 

Tile kinetic theory of neutral gases, as it is presently 
understood, is based on the work of Maxwell and Boltzmann. 

The assumptions made by all previous workers, namely, that all 
molecules of a gas move with the same speed, was abar.doned 
and the random character of molecular motion was recognized. 

The Maxwellian distribution of velocities for a uniform gas 
in equilibrium, and the law of equipartition of the mean 
particle energy in a mixture of gases were firmly established 
[2,7]. In an attempt to give a rigorous justification of 
Maxwell's assumptions, Boltzmann derived an integro-diff erential 
equation to describe the evolution of the velocity distribution 
in space and time. This integro-diff erential equation is nov; 
known as Boltzmann equation. 

2«lol The Velocitu-Distribution Punction 

3 

A gas particle can be characterized by its position r 

Jk 

in space and its velocity v or momenta p. ¥e will be using 
V or p according to convenience. Information is not only 
required about the spatial distribution of particles, but about 
their velocity and momentum distribution as well. Only with 
the latter we can study the momentum and energj'' flows. 
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The simplest function that contains the desired infor- 
mation is the velocity distribution function, f(r, v, t) , 
defined in such a way that f(r, v, t) drdv is the expected 
number of particles, at time t, in a volume element dr located 
at r, whose velocities lie in dv about the velocity v. The 
volume elements dr and dv are, by defiii tion. 


dr 

5 

= d^r = dxdydz, 

(2.1) 

dv 

= d^v s dv dv dv 

X y 2 

(2.2) 


where x, y and z represent the cartesian coordinates and 
v^, Vy and v^ are components of v in directions of x, y and z, 
respectively. The macroscopic or average properties of the 
gas as well as the flux vectors of mass, mcmentum and energy, 
can be defined in terms of the velocity distribution fimction 
(see Chapter 3 ) « 

2olo2 The Boltzmann Equation 


Boltzmann derived the kinetic equation governing the 
evolution of the velocity distribution fimction in space and 
time. Many modern derivations have been developed by several 

„ JJ|. jj, 

authors [2, 7, 1]. ' equation for f(r, v, t) can be 

written in the form 


df _ / 3f \ 

dt ~ ^ a t ^ c 


( 2 . 3 ) 



10 


The left 

hand side of 

Eq # 

2.3 

is given by 


df 

dt 

Jf 

3 t 

dt 

. if 
dr 

dv 

+ — 

dt 

. il 

-i 

3 V 



3f 
" 3 1 

jk 

+ V • 

V f 
r 

m 

• V f, 

(2.4) 

wh er e K 

is the 

external force. 

m is the mass 

of th e 

individual gas particles. 

and 



3 

3r 

= V 

r 

= i 

Fx 

3y 

+ k ^ , 

3 Z 

( 2 . 5 ) 

and 

3 

9 V 

- % 

= i 

3 

,*.3 t 3 

(2.6) 


4> A. A 

where i, j and k are the mit cartesian vectors. 

The explicit form of the Boltzmann collision integral (" 3 ^)^ 
is not of interest to us here, since the assimiptions mads in 
its derivation are not applicable to plasmas as discussed in 
ChapteT' 1. We will, therefore, proceed to a kinetic descrip- 
tion o± plasmas in the next section. 

2o 2 The BBGKY Hierarchu of Equations 

Most of the kinetic theory of plasmas is based on one 
of the following three approaches : 

(a) Pokker-Blanck approach, 

(b) BBGKI hierarchy of equations, and 

(c) The method of diagrams. 
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In the Pokker-Planck approach, the statistical method 
for ti-'eating random processes is applied to ions or electrons. 
The fact that the prohahility of a large engle deflection of a 
charged particle in a plasma oocuring as a succession of many 
small angle scatterings is much more than as a sin^e binary 
scattering is used as the basis for the Pokker-Planck theory. 
The starting point in the other two approaches is the liouville 
equation . In the method of diagrams, developed by Prigogine 
and Balescu [24,25], the Liouville equation is Pourier 
analysed and sol-zed by iteration. In the EEGEJ (Bogoliubov, 
Born, Green, Kirkwood and Tr on) method, a hierarchy of 
equations is derived from liouville equation, which are then 
solved by successive approximations [2j. Since the work 
reported here is based on this method, we will discuss BBGKI 
system of equations in scaie detail. 

The iouville Equation 


The microscopic state of a plasma or any other many -body 
system is completely defined by the coordinates and velocities 
of the particles. Similarly the statistical properties of 
such a system are completely determined by the spatial and 
velocity distribiition of its particles. If a system has N- 
particles, the K-T)article distribution function, is defined 

as follows: 

P» ^2 **•' ^1* P2****^K* * 


(2,7) 
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where q and p denote the coordinates and the momenta of 
the N-particles, and is such that 


V ^ P,j(q,p,t) 


'N 


dqdp = 


V q^p y t ) dq-j^o.q2* • • dq-j^^dp^^dp 2* • • dp|^) 


A 


3. ^ Jk 


( 2 . 8 ) 


is the probability that at time t the coordinates and 
velocities of the particles have the values q^, qp***qjf» 
Pl, P 2 ****Pj;j ill iiie range dq^, dq 2 ...dqj.-, dp^, dp 2 . ••dp^j., 
respectively. is normalized according to the following 

relation % 


. A 






rB 


(2.9) 


The evolution of P^^ in time and phase-space is governed 
by the liouville equation [4,2]; 


dP 

Tt 


N 


3P 




+ 


H 

I ( 

i=l 


sq-i '^i ^ 


- T*Pi) 

1 


0 

(2.10) 


where a dot on the top of the variable denotes derivative 
with respect to time. 

By means of Canonical equations of motion [23]» 


^i = 




s 

P 

m 


Pi 




i = 1,2,. . . ,B', 

aq^ 


( 2 . 11 ) 
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where stands for the Hamiltonian of the H-particles, 




= I [Tj- (Kinetic energy) + V^(Potential energy)] 
i=l 


( 2 . 12 ) 


We can write the Liouville equation in the form 


ap 


IT 


9t “ ^ %* % 


(2.13) 


or, 


3P 




at ~ 


(2.14) 


where IT 

^ SH.. ap 3H„ 3P„ 

[%,%] = I - -f~^) 


(2.15) 


is the P oisson-hrackot expression for and P^^. The dots 
between the derivatives on the r.h.s of Eq. (2.15) denote 
a scalar product of two vectors. 

Por an external force-free condition, which will he 
considered in this work. 


E 




N 2 

I I 

i=l 1 < i < 3 < E 




( 2 . 16 ) 
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where | Pj_t ? and 9'^^^ is the interparticle potential, 

a function of interparticle distance and not of the individual 
coordinates, that is, 

(2.17) 


¥e define the following operators for an external 
force-free system of H-particles ; 

N E 


% = J I 

^ 3q.i 

i=l i=l 


V . 
1 


V , 
1-: 


(2.18) 


E 


Q = y (_iii 




-i 


(2.19) 


3q. 8p. 


Sqj SPj 


1 _< i < 3 <_ h 


1 

2 


E 

I 


E 

r ^ 0- ■ a 


i=l 3=1 




ap- 


aP- 


( 2 . 20 ) 


In terms of these operators the liouville equation can he 
written as 


IV + ( - V = 0 (2-2« 

where (K^-Ojj) can he identified as the Hemiltonian operator 
for an external force free-system of N-particles [7 ]. 
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The BBGKY Equations 

For brevity of notation, we define p^) , 

i In order to obtain the BBG-KT equations from 

the Liouville equation, we define the reduced distribution 
functions, Fg(x^. . . x^jt), s = by 


i-gUi.Xj. 


Xg,t) = 


V 


-(N-s) 




"S 


,t)dx. 


!+l* 


.dx, 


•N 


( 2 . 22 ) 


where V is the total volume of the system. In order to 
concentrate on volume effects and minimize boundary effects, 
the system is assumed to be infinite by allo\d.ng both II and V 
to become infinite in such a way that the particle density, 
n = E/V = finite. 


By integrating the Liotjville equation over 
it is strai^t forward to obtain the so-called BBGEI hierarchy of 
of U equations, which were derived independently by Bogoliubov, 
Bom, Green, Kirkwood and Yvon [3»5,6]. 


9F 

s _ 

9t “ 



I'.c,] +^/dx,^T [ I 0, 

i=l 


g. ■ V •» ^s+1 ^ ^i,s+l’"- s+1- 


(2.23) 


s = 1, 2, .... ,11, 


where the square brackets denote the Poisson brackets (see 

Bq. 2.15). For a derivation of Bq. (2.23), reference should be 

made either to the originM.- works cited above, or to Ref s. [ 7 26]. 
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Again considering external-force free system and by 
defining the operators s 

V Pi 3 

^ ^ - (2-24) 

i=l ^ s'li 


and 


e 


1 1. i ^ j £ s 


(__L _ \ ) 

9q.i 3 Pi 3Pj 


(2.25) 


e, = 0, 


(2.25a) 


s 


Z / dx . — j— i*-— , ;j=s+l: 

8 li aPi 
i=l ^ ^ 


(2. 26) 


and taking the limit of (H'-s)/V' = n as N and Y tend to be 
infinite, we get the following alternative form of HB GEf 
equations ; 

8P 

TV + , (2-2U 

B=l, 2 , . . ,H. 

As pointed earlier (K - Q ) can be identified as the 
Hamiltonian operator for the external-f orce-free system of 
s particles and is sometimes called the phase mixing 
operator [7 ]. 
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Ths Generalized Kinetic Equation, Irrevercibilitu and 
Statist ical Independence 


The equations (2.27) are equivalent to the Liouville 
equaoion, since they have been obtained from the latter 
withont making any additional assimnptions. Hovrever, it can 
be seen immediately that equations (2.27) form a chain : for 
every s, the equation for contains the ■unknown The 

first of these equations, namely 


Pi 


^3t m-i 
1 


9 0, 

n / / dq^dp 


3 F 

12 ^ 2 


3qi 3Pi 

( 2. 28) 


can be called a generalized kinetic equation, applicable to 
any external-f orce-free system. The left rumnber of this 
equation has exactly the same form as the left hand side of 
the Boltanann equation. On the other hand, the right hand 
member of Eq. (2.28) gives an exact expression for the 
change of due to interactions. It cannot, ho-wever, be 
solved without making additional assumptions, because ^2 
unkno-wn. If we write the equation for Fp, it will contain 
the unknown F,,, and so on. 

Another very important feature of the BBGKJ Eqs. (2. 27) 
or the generalized kinetic Bq, (2. 28) is that they are never- 
sible in time, just as the Liouville equation from which they 
have been derived. Therefore, unless a definite direction of 
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time is introduced into these equations, they will not describe 
the irreversible approach to equilibrium from non- equilibrium 
states. 

It may also be noted here that in terms of the reduced 
distribution functions, statistical independence of the particles 
can be defined as follows. If 

s 

I’g(x^,...,Xg,t) - n F3 _(x., t) = 0, (2.29) 

i=l 

s ^ 2, 

the particles in the system are statistically independent 
or uncorrelated. This is also called the assumption of 
molecular Chaos [16], This is one of the assumptions made in 
the derivation of Bolt 22 nann equation. Since in a plasma, due 
to the long range nature of the Coulomb force, the particles 
are correlated i. e. Eq. (2.29) is not satisfied, a different 
approach is called for. 

In Section 2.5, we present briefly, Bogoliubov's method 
for solving the BBGhT equations, our main interest being in 
the equilibrium one-particle and two-particle distribution 
f -unctions for an external-f orce-free fully ionized plasma. 

In the first order approximation Bogoliubov's method solves 
the BBGEy equations for small values of the plasma parameter, 
e, which is defined in Section 2,4* 
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Z&g— fi.CSA^__.ParaPterer and the Plasma Aoprozcination 


Tile plasma parameters e, is a measure of the inverse 
of plasma particles in a Debye sphere and can be defined by 
[14] . 


e 


1 


n r-r 


(2.30 


1 


where K is the particle density, and rj^ is the Debye length 
given by [26, 27, 14, 2] 


e kT 1/2 

( -V ) 


D - '2 

ne z 


e © 1/2 

o \ ' 

2 ' 
ne s 


(2.31) 


In the last equation k is the Boltzmann constant, e is the 
of 

absolute value electron charge, T or 9 is the temperature 
defined in terms of the velocity distribution function in 
Chaptex 3, and is the permittivity of free space 
(Rationalized M.K.S, or SI units are used in this work). 

It can be easily seen that the plasma parameter, s, 
is also a measure of the ratio of the mean interparticle 
potential energy to the mean kinetic energy of the particles 
in a plasma, that is 


< PB > 


e -v 


< iCB > 


( 2 . 32 ) 
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If e 4.^ 1, the interparticle potential energy will be much 
less than the mean kinetic energy of the particles. This wU-1 
be so far a sufficiently dilute plasma. The assumption c ± 
is called the plaana approximation. It may be noted that as 
the plasma particle density decreases, the number of particles 
in the Debye sphere increases, and then e decreases. 

5 Boooliubov's Approach to the Solution of BDGKY Equations 

The BBGKI hierarchy of equations can be solved by three 
different approaches s 

1 . B ogol iub ov ’ s M eth od , 

2. Frieman and Sandri Method, and 

3. Cluster Expansion. 

As was mentioned in Chapter 1, the three approaches are 
equivalent. All the methods, in their first order approxi- 
mations, solve the BBGEJ equations for small values of the 
plasma parameter. We will discuss the Bogoliubov*s method in 
some detail. 

The assiimptions made in Bogoliubov’s theory are ; 

1. The evolution in time of all s-particle distribution 
functions, E , s ^ 2, can be described through their functional 

dependence on the one-particle distribution function 

* * * ^^s*^l^ ^ 


Eg(x^,. .. ,Xg,t) 


s ^ 2, 


(2.33) 
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With the help of this assumption the generalized kinetic 
equation can be written in the form 

» (2.34) 

where A(x^ is a functional of 

2. The past is defined as that directi on "of time in vrhich 
correlations among the particles decrease to zero : 


[Fg(x^. . .Xg,t)- 


n F^(x^,t)] = 0, 
i=l 


(2.35) 


s > 2. 


3. The following expansions in terms of the plasma 
parameter, e, are made : 

F„ = I s > 2. (2.36) 

s s 

rt^o 


I ( 2 . 37 ) 

n=o 

When we utilize Bqs, (2,35) to (2.37) in the BBG-ET equations, 
we obtain equations to each order of e, which can be solved 
in succession. 
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2o 6 The Eauilibriun One-Particle and Two-Particle Distribution 
Tunc t i ons 


V/hether BogolialDOv ' s theory or other equivalent 
approaches describe adequately the irreversible approach of 
non- equilibrium plasmas (or other systsns) to a uniquely 
defined equilibrium condition is not yet kno^m. Hov^ever, to 
the first order in e, it is known to satisfy the Boltzmann 
H-theoren [7]. As an extension of this fact, it is generally 
assumed that the Bogoliubov’s theory does describe an irrever^ 
sible approach to equilibrium. Guernsey and len£'rd[ 2] have 
obtained the following one-particle and two-particle equilibr: 
distribution functioiis, to the first order in e ,2 


Fn(x,,t) = P®(p) = (2ii:mG)“^/^ exp(-p^/2me) 

Equilibrium 

(2.58) 


and 

E.(x. ,X5,t) = F®(r,p.,P2) 

^ ^ Equilibrium 

e2 

= 5'^(P]_)^' ” 4n e ^r0 ®^P^” ^ ^ 

(2.39) 


where 

m = mass of an individual particle, 

0 = temperature at equilibrium, 

^1® ~ charge of particle 1, 
z^e = charge of particle 2, and 

r = ^ distance between particles 1 and 2 

= fVW* 


vm 
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Z-| Z p6 

The factor [l - exp(- --) ] in Eq. (2.39) shows the 


47teQre 

effect of interparticle correlations in a plasma. V/e will 
discuss it further in Chapter 4. It should be noted that as 
required by basic definitions F®(p) in Eq. (2.38) normalizes 
to imity when integrated over p (see Appendix I). 


¥e now proceed to Chapter 3? where we define the 
macroscopic variables, the flux vectors, the macroscopic 
conservation equations, the rate equations, and the transport 
coefficiaits. Then in Chapter 4, we discuss the hydrodynamic 
approximation for the distribution functions and obtain 
appropriate expressions for the transport coefficients 
considering small deviation from equilibrium. 



Chapter - 3 


llACROSCOPIC CONSERVATION E >UATIO!JS AND 
TRA NS PORT COEFi 'EC I ENTS 


3^1 The Eacroscopic Variables 


Tile macroscopic variables of interest in this work 

^ •A . 

are the particle density n{q,t) or the mass density p (qjt), 

JL At 

the macroscopic or mean velocity uCqjt) and the temperature 
TCq^'fc) 01’ ©(qjt). ¥e wish to define these variables both in 
terms of the velocity distribution function f(q, v, t) and 
the reduced one— particle distribution function p» tj 

introduced in Chapter 2. 

The following relationship between fCq, 'v , t) and 
p, t) holds from their respective definitions (Sections 
2.1,1 and 2.2.2); 


f(q,^.-fc)dT 

= n S' 3 _(q,p,'t)dp , 

(5.1) 

where n = N/V, the 

particle density- 

- 

-i , 

The macroscopic variables n, ^ ,u and 

T are then 

.a -i 

straight forwardly defined in terms of fCq» i t 

t) as 

tl 

f (q, v,t) dij- , 

(3.2) 

p(q,t) = mn(q,tj ^ , 

(3.3) 


= -i— / ^ 

nCqj'fc) 


(5-4) 
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and 




1 t 1 

— j ^ 

n(q,t) 2 


mU^ 


f (q,v,t)dv 

(3.5) 


where U = jU! = }v-u | . U is called the thermal or peculiar 
velocity. Extension of Eqs. (5.1) to (3.5) to systems of more 
than one species is straight forward. It is easily seen that 
the equilibrium distribution function of Eq. (2.38) satisfies 
above definitions with u = 0 (see Appendix I). 


Jo 2 Flux Vectors 


The flux vectors of interest here are the particle 
flux, r (q,t), the momentum flux (which is actually a tensor 
called pressure or stress tensor) p(q, t), and the heat flux, 

J, defined as follows i 

r(q,'fc) = / ^ f(q»v,t) dv 

= n(it) u(it) (3.6) 

f = [PiJ i,3 = 1.2.3, 

= m / f(q,v,t) dv (3.7) 

j = |m / U2 U f(t^,t) dv (3.8) 

All the quantities appearing in EqS. (3.6) to (3.8) have already 
been deiined. It follows from the defination of p and T that 
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P^i = nkT = ne, i= 1,2,3. (3.7a) 

It should bo pointed here tha-t in a system with 
interparticle potential energy, the pressure tensor and the 
heat flux will, in gc-neral, have terms relate)! to inter- 
particle potentials. These terms have hetai derived by Chch 
and Uhlenbeck [2 ] and are very complicated to be quoted 
here. In the plasma approximation, e < < 1, these terms are 
negligible. 

Th^Ilacroscopic Conservation Eouationa 

The macroscopic conservation equations are obtained 
by multiplying the kinetic equation by the conserved quantities, 

wit 

namely mass m, momentum mU and the kinetic energy -5-mv^, in 
turn, and integrating over v. If the potential energy flow is 
iiegligible, these conservation equations talce the following 
form for an external-f orce-free system [2,16] 

■^ + v-( r ) = 0, 

or , 8 p 

^ + V -(mr ) = 0 , (3.9) 

P + U* V + V .p = 0, (3.10) 

nd _ 

[ ^ ^ + fev -J + pJ] = 0, (3.11) 
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where D 


is the strain-rate-tensor , whose component are 



n SU. 3U. 

:i( i. . — J 


(5.11a) 


The sin^e dot product between a vector A an<i tensor B 
is a vector whose components are 


-i _ 


^ ^3 ®3i* 


The double dot product between two tensors A and B is a 
scalar given by 


B = 


3 


^i3 ®3i* 


of 

The extension^ Eqs. (3. 9) to (3.11) to systems with more than 
one species is strai^t forward, provided the flux vectors 
defined in Section 5.2 are also rewritten for individual 


sp ecies. 


3 o Th e Plate Equations and the Transport Coefficients 

The phenomenological rate equations for transfer of 
heat, momenttim and mass (or number of particles) are : 

1. Fourier’s Equation of Heat Conduction s 


J 


■X 7 rp 


( 5 . 12 ) 
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2. Newton's Equation of Viscosity % 

Pij = - ^ . i+3. (3.13) 

3. Eick's Equation of Diffusion t 

r = “ D V n. (3.14) 

Equations (3.12) to (3.14) define- the transport coefficients 
of interest in this work, namely, the coefficiait of thermal 
conductivity, Z, the coefficient of viscosity, ji., and the 
diffusivity or diffusion constant, D. It should he emphasized 
that we have assumed Z, and D to he independent of direction. 
This will he the case if plasma is free from external fields. 

In Chapter 4, the hydrodynamic approximation of the 
distribution functions and the generalized kinetic equation will 
De considered. It will he seen that if we retain only the 
first power of the gradients of the macroscopic variables, the 
moments of the kinetic equation yield the rate equations (3*12) 
to (3.14) together with the appropriate expressions for the 
coefficients Z, \x and D, defined above. The expressions obtained 
are fomd to be in good agreement with the results of other 
authors based on more advanced theory, for example E okker-Planck 
approach [21,22] or a^ p erturbati on solution of BBGZI equations 
[ 2 , 20 ]. 



Chapter - 4 ' 


TilA^^UPuFT COEFI' I-Zini’TS IN THE HYDBODYrAhlC aPFRO.:Ii' 


■.TIOR 


T.hF._liiidrodun.ani c Y^DDr'ozzi.ia t 


ion 


In the hydrodynamic approximation [2,20] to the kinetic 
equation, it is assumed that the distribution functicnsdo not 
depend on space and time explicitly- The spatial and time 
dependence is assumed to he only through a functional dependence 
of the distribution functions on macroscopic quantities n(q, t), 

•Jk «JI 

u(q,t) and T(q,t). That is, we assume for (= F) and I ' 2 'the 
following formss 

F(x,t) P(S I n, u, T) (4,1) 

and 

^ ’ ^ 2 .* ^ 2 ^ ^ 2 ^ ^ 2 ^ * 

(4.2) 

where rp2 = Ki--q2^ 

It should be emphasized here that the functional depen- 
dence on n, u and T is understood to mean a dependence on n, 
u and T as well as on any of their spatial derivatives [7]. 

To state it explicitly, let p be a 5 -cQmponent vector including 
n, u and T. The F and F^ can be written as 

F(x,t) F(p, p, ^ p, )' (4.3) 
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and 


^ 2 ’ ^ ^ 2 ^ ^ 12 * ^ 1 ^^ 2 ’ ^^ 2 ^***^ 


(4.4) 


where the dots indicate hi^er order spatial derivatives of p. 
The iterative solution of Chapin an-Enskog [l] and Choh-Uhlenheck 
[2] for neutral gases and of Sund art^an-Wu [20] for fully 
ionized plasmas are b aued on the assumptions ©nbodied in Bqs. 
(4.3) and (4.4). 


4olol Limitations of the Hydrodynamic Approximation 


The hydrodynamic approximation is valid only if the 
system is not far from local thermodynamic equilibrium. 
Quantitatively, it means that the length over which fractional 
change in the macroscopic variables is sigrificant is much 
larger than the mean free path, X of the particles, that is 


X >1^ 

n 


< < 


1 , 


(4.5) 


9Ui ^ 

X — i-i< < 


(4.6) 


and V T d < i 

X T 


(4.7) 


Under these conditions, a particle will encounter approximately 
the same values of the macroscopic variables in its one free 
flight. The hydrodynamic approximation can then be mathematically 
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justified based on Hilbert's uniqueness theorem [7 ]. ThuSs 
for a system not far frcm local thermodynamic equilibrium, the 
ini tial-valu e problem posed by the kinetic equation can be 
effectively replaced by an initial -value problem in hydrodynamic 
equations. 

It should be pointed out here that retaining hi^er™ 
order derivatives in Eqs. (4.3) and (4.4) or hi^er powers of 
first order derivatives, does not make the solutions obtained 
by this method applicable to systems with large gradients. It 
only refines the solutions further. More moments of the kinetic 
equation should be included if the system is far from equilibrlur. 
(see, for example Grad's method of thirteen moments [2]), 

The above brief discussion provides seme insist into 
the applicability ox the concept of transport coefficients. 

Even if non-linear rate equations are used to define these 
coefficients, valid solutions can be obtained only for plasmas 
not far from local thermodynamic equilibrium. To illustrate it, 
we can cite the example of laser-pellet-fusion. The plasma 
created on the pellet- surface is far from local thermodynamic 
equilibrium. Th'as in our view, no valid hydrodynamic model based 
on the five moments n, u and T can be developed. Either the 
kinetic equation should be solved as an initial value problem, 
or more moments should be included. 
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^ <>2 Characteristic tlean Free Path for a Plasma 

In a plasma, due to long range of Coulomb force, there 
is no well defined mean free path, because particles interact 
all the time. However, to define a state not far from local 
equilibrium throufeu Eqs. (4.5) to (4,7), we can choose a suitable 
characteristic path length, for example, a path length over 
which a particle, by random walk, scg-.ttors throu^ a total an^ e 
of 90^ by multiple collisions. This path length for electrons 
in a fully ionized plasma is given by [14 j s 

(e^/^e^kT)^ fnA]“^, (4.8) 


where Ze is the charge on ions and 


A = 


3 V2 

2 '■ Ttn ■' 





0)V2. 


( 4 . 9 ) 


(|n A is called the Coulomb logaritlim. 

Over a large range of kT and n, iin A is found to be a weak 
function of plasma parameters [14]* 3)he characteristic path 
length for ions, is 'obtained upon I’eplacing Z by Z in 

Eq. (4.8). In writing Eq. (4.8), the Co-ulomb logarithm is 
assumed to be same for electron-electron, ion-ion and ^ectron- 
ion collisions. This is a reasonable approximation due to slow 
variation of ^n A given by Bq. (4»9) actually applies 

for electron-ion collisions. If this assumption is not made, 


V 
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plasma ion and electron densities will appear separately iis. 
Eq. (4.8), and the corresponding expression for and 
can he easily written. 


4 o .5 Non-Qaui librium Distribution Functions 


(a) One-Particle Distribi^tion Punction i 


In accordance with the hydrodynamic approximation, xf 
the plasma is not far from equilibrium, the spatial and temporal 


variations of non-equilibrivm distribution functions are 
adequately described by their dependence on 0 , i\ and n, and on 
their spatial gradients. For further development, we make the 
following assumptions s 


nrn X ,'\yn 

1 . Since the quantities A 7 ^ and A — ^ are 

i ^3 

small ( 1 ), their higher powers can be neglected. 


2. Without loss of generality, the direction of a 
gradient is taken to coincide with the x— direcuion. 

3 . Again, -vjithout loss of generality, the gradients in 
macroscopic variables are treated one by one to study 
the particular transport process associated with that 
gradient. For example, to study thermal conductivity, 
a temperature gradient is assumed to exist, while 
gradients in the macroscopic velocity or density are 


taken to be zero. 
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With the help of these assumptionsy the one-particle 
distribution function for a non-equilibriun situation can be 
written as 

= ?^(p)[l “ (4.10) 

where F (p) is the Maxwellian distribution given by Eq. (2.38)> 

p is to be identified with any component of p , in particular, 

6, u and n respectively, and h- (p) is to be obtained by 
y p 

substituting Eq. (4.10) in the kinetic equation. The negative 

sign of the last term in Eq. (4.10) has teen choosen for 
convenience. The function hp(p) must satisfy the following 
properties ; 

1. It must be a dimensionless function of p, i . e. 

hp(p) = hp(i) (4.11) 

where p' = p/Vm© (4.12) 

2. It follows from basic definitions that 

00 

/ E-dp = 1 

- 00 

Since 

/ E® dp = 1 

00 

I (see Appendix I ) , we must have, 

CO 

/ h.e(p) <ip = 

- CO H 


0 , 


(4.13) 
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Obviously, iHq. (4.13) is also satisfied if we repls-ce p by p'. 
(b) Two-Particle DistribiAtion Pune ti on; 

It can be easily seen that the generalized kinetic 
equation applicable to plasma (Bq. 2.28) contains not only P^, 
but also the two~particle distribution function P 2 " Hoither 
the solution nor the moments of the kinetic equation can be 
obtained without knowledge of P^. By analogy with the equili- 
brium two-particle distribution function (Bq. 2«39), we assume 
that for a situation not far from equilibrium, ?2 can be 
written as 


Jt — ' — 1 V - 

2^Pq»P2* 2° * ’ * ^ ~ ^(Pq>Pq: 


_ .i _ — 

)! (p 2?^ 2 * * * * ^ Gr(rs0j n) 


(4.14) 


where P is given by Eq. (4.10) and 


G(r,©,n) = 1 


Zj_Z2G 

4^ e ^r© 


sxp (■ 


(4.15) 


r- 


D 


where, 

^ = "*12 
© = ( e ^+© 2)/2 

n = (n^+n 2 ) /2 . 

r^j is to be obtained by replacing © and n in Eq. (2.31) by © 


and n. 
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It may be pointed out here that “Pg ^ singularity at 
r = 0. This can not be removed unless quantum effects are 
accounted for. P or a detailed discussion reference should be 
made to [2]. It should also be emphasized that the equilibrium 
form of F 2 (Eq» 2.39) is applicable only for dilute plasma 
(e < < 1). Hence the results obtained by taking small deviations 
around this form will also be applicable only to dilute plasma. 


I 


Even with the above simplifying assumptions, it can be 
easily seen that the integral on R.H.S. of Sq, (2.28) can not be 
performed analytically. If V©(= kvi) and 7 n are small, the 
function G(r,©,n) can be expanded in a Taylor series around 0^ 
and Again only the first order terrr-s in the gradients are 

retained. The detailed expansion is given in Appendix II. 


^ Collision Term in the Generalized Ki/tcti c Equation 


With the approximations made about the form of 2 ^ 

non- equilibrium situation in Section 4.3b nd Appendix H , the 
collision integral in the generalized kinetic equation, Eq. (2.28), 
can be easily evaluated. Considering 03aly a fully ionized two 
species plasma, we have 



^ 1 ^ 2 ^ 

4Txe^r 


(4.16) 


where is the charge on particle 1 and Z 2 e is the charge 

on particle 2. It should be noted that the particles 1 and 2 
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could iDelong either to the same species or different species. 
Thus there will be two terms in the collision integral? ootained 
by fixing Zj^ and varying Z 2 over the two species present. 

This summation ca.n be performed at the end, after the collision 
term has been evaluated. 

Denoting the collision integral in Eq. (2.28) by I and 
utilizing Eqs. (4.14) to (4.16), we got, 


90 


n 


I dq G(r,e,n) / dp 2 5'(2) 


?Pl 


'll 


(4.17) 


where F(l) and F(2) denote F at points 1 and 2 respectively. 
By definition, the last integral in Eq. (4.17) is unity. Tne 
other integral in Eq. (4.17) is evaluated in Appendix III. 

The result is 


Tca 


■■D 


2 2 


30 


V 9 
© 


+ 


y n 


)] 


n (4.18) 


2 

where a = e 

The subscript 1 has been dropped, since it is no longer necessary 
The generalized kinetic equation for electrons now takes the 
following form s 


3F 
9 t 


+ 


P ^ 


m 


e 3q 


2 

na r- 


D 


3e 


(n +11. Z^) ^ 
' ® ^ 3p 


-V© 

"G 


TIB. T-qIX Z 

______ 


9 p 


fV© + vn -1 

'l-©~ ^ 


7 n *! 
n 


(4.19) 
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■where is the electron-mass, n and n- are the average 

c* -i- 

electron and ion densities, respectively, and ali othen quantities 
have already been defined. In the last £•', ep above, we have used 
the equality 


n 


+ nj_ Z 


= n Z 


(4.20) 


which is readily identified. 

It can be easily seen that the icinctic equation for 
ions can be obta,ined by replacing m ^ by m. and Z by in 

O 1. 

Eq. (4-19). 

By substitiiting F frem Eq. (4-10) into Eq^ (4.19), the 
governing equation for h^ can be obtained. It should be noted 
that with F from Eq. (4.10) 


3 F ^ 5 B 

3 t - 3 p *3 1 ’ 


( 4 . 21 ) 


iZ _ nr 

^ - 3 p * ^ 

3 q ^ 


( 4 . 22 ) 


Since derivatives of p with respect to time or space can 
be neglected as higher-order terms. The resulting equations 
for hp are then to be solved together -s-ri-th the macroscopic 
conservation equations and other defining equations for macro- 
scopic variables and flux vectors given in Chapter 3. It is 
seen in the next section that h^ appears in the transport coeffi- 
cients only through a definite integral. Since h^ is a 


dimensionless fimotion, this integral efiects the transport 
coefficients only by a multiplicative constant, which will 
not be evaluated in this work, 

f of Expressions of Transpor t Cog fric ioj} t s 

The expressions for the transport coefficients can now 
be obtained by taking mem ents of Eq. (4ol9)o This means multi- 
plying Bq. (4.19) by the particle density'' (n or n. ) , mcmentum 
n^mj^U) and the kinetic energy or , 

in turn, and integrating over all p. The following identities[ 2] 
are utilized in obtaining the manents i 


jt 


/ Q(S) 

P 3? 

= 

^ ’/ Q(p) 1 S' 

(4.23) 


tv iw- 

2 q 


/ Q(p) 

dp = 

- / dp, 

(4,24) 


^ p 




where Q(p) is any property associated with p. In writing 

Bq. (4.24) it is assumed that Bo -> Oasp-^+®, 

4o5.i Ion Hiffusiuity 

Even thou^ electrons are more mobile than ions, it is 

the diffusion coefficient of ions that governs net diffusion 

in plasma. This is due to ambipolar effects [27]. To investigat 

the ambipolar effects, the kinetic equation has to be solved in 

conjunction with Maxwell’s equations. It can bo proved that if 

T - Tj and m. > > m , ambipolar diffusion coefficimt in twice 
Q i 1 B 

the ion diffuaivity [27 ], 
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To obtain an expression for ion diffusivity, we write 
Bq. (4.19) for ions with 


= 0, vn = 0, VTx iz 0. 

Under these conditions, the kinetic eoustion is 


mi 

9 p 

Jbe 

9q 

A 

iJl 4. rL. . 

9n ''it ^m^ 

"p’ rn 



-ita^rj^n ^ 

“ 3© 

9p 

711 

n 


(4.25) 


Multiplying Eq. (4.25) by Q = n^, and utilising Eqs. (4.23) 
and (4.24), we get Eq, (3*9) for the conseinration of particles. 
The x~cQmponent of the particle current r is given by 
Eqs. (3.6), (3.1) and (4.10) as 

^x ^ ^i vn)dp 


•Hi X . gn. 

= /“• ^ Php)[i hjp)]ds 




n^ ^x 


n' 


X. 9n. “3/2 00 

— — (atm^e) I p^e 

3 - 00 


*p^/2m^0 


h-J^(p)dp 


( 4 . 26 ) 


^ -A 

since f p E (p)dp = 0 (see Appendix I) 
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By changing the variable of integration from p to p* = p / y'inS' 
and noting 

a 

dp = (m^e) dp' , (4.27) 


and using 


A ^ -LX VJiil 






r = ^ 
z 


9e 


0 




5/2 


(Z+1)Z^ £nA ][m^ n^^ ^ 


(4.28) 


where we have used the fact that n = n^(Z+l)j and is given 


I 


n 


CO 



-(p')^/2 



(4.29) 


Comparing Eq. 



(4.2.3) >7ith the rate Eq. (3.14) 

In ( &5/2 , 

V"'^(Z+1) Z^ A 


we get 

(4.30) 


Since is a dimensionless quantity, the functional dependence 
of on plasma parameters is given by the last term in 
Eq. (4.30). 

It can ij^ GGon from Bq.(4.30) that the diffusion 
cocficient is ixlv jrsely proportional to the nquara root of mass. 
Hence tK«,> electron diffusivity would be much larger thtJl ion 
diffusivitysHowever, because of ambipolar effects, the much 
slower ion diffusion governs the net diffusion rate in the plasma. 
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^ho5->2 Ion Viscosity 

At or near equilibrium in a fully ionized plasma, the 
ions and electrons carry approximately the same amount of energy 
per particle. However, since the ions are much more massive 
than electrons, most of the momentum is carried by ions. Hence 
viscous effects are mainly due to ions. In this section, we 
will develop an expression for ion viscosity. 

For this purpose we assume that the density and tempera- 
ture are uniform throughout the plasma, and that the deviation 
from equilibrium is due to small gradients in the macroscopic 
velocity , i. c . , 


V n = 0 = V ©, 


Vu 4 0. 


Further, to choose a one-dimensional situation, we assume vu^ 4 0, 


V u 


9 

0 




Under these conditions, the kinetic equation. 


Eq. (4.19), along with Equations (4.21) and (4.22), becomes 


9F. 
3 t 


JUi 

P 

“i 


IL 

3 1 


3U 


^x 3F 


3U 

Ljl + 

y “l 


a u 

3^ 


= 0 (4.3l) 


It is straight forward to sec that by multiplying Eq. (4.31) 
by Q = n^m^v, and integrating over all p, wc obtain Eq. (3.10) 
for the conservation of momentum. Iho components of the stress 
tensor aro given by Eq. (3. <) » 
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P 


xy 


1 


m. 




( 4 . 32 ) 


where 


^ = 


Jjju 

nij^U 


m . 
1 



(4.33) 


utilizing the assmued form of P for near—ea'oilihrium conditions 
from Eq. (4,10), we obtain, 


xy 


^i /• T? ^ "i ^ A ^ 

m . / ^ ” u” 9 X 

^ y 


(4.34) 


It should he noted here that in the presence of macroscopic 

P 

velocity, P is to he taken as ; 

F® = (2icme)"^/^ exp(- ^^/2me) . (4.35) 


It is then easily seen that 

CO 

/ ^ P^ dp = 0 (see Appendix I) 

—■00 ^ J 


Thus we have. 


■D 

' xy 


n . 


m. 


1 1 

u. 


y 


9U^ -5/2 

(pjxmj^©) 


/ >?C; ^y 

CO 


/2mj_e 




( 4 . 36 ) 
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By defining the nondimensional quantities, 

j. ^ 

>' = 7^ • "y = 

and noting that 


It 

(m^e)^/^ 


We obtain, 




n . ^ . 

Vm .0 I 

P 

1 a 

* X u 

xy 



(4.36ft) 


where I is the diniensionless integral given hy 




u u. 


(4.37) V 


y r»» 


Using from Section (4.2)* and noting that n = nj_(Z4-l), 
we get 


__9< \ iTm. 

^^y Z^(Z+l) ^nA 


(4.38) 


Comparing Eq. (4.^) with the rate Eq. (3.13), we obtain the 
desired expression for the coefficient of ion viscosity: 

9e? I fm, 


p = 


Z^(Z+l)enA '"e^ 


(4. 39) 


It can be easily seen that the viscosity is directly proportional 
to the Square root of mass* Hence viscous effects due to 
electrons will be negligible. The expression derived in Eq. (4. 39) 
is in good agreement with that given by Eokker-Planch theorjCSf 21]. 
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4o5»3 Thermal Conduct ioitxj 

The phenomenon of thermal conduction in a fully ionized 
plasma arises mainly due to collisions of electrons x/ith ions 
or other electrons. Hence we will develop in this section an 
expression for electron thermal conductivity. The method is 
very similar to that used in Section (4.5.1) and (4.5.2) for 
diffusivity and viscosity. Here, v^e assume that the departure fr 
from equilibrium is due to a small gradient in t emperatur e alone, 
i. e. , 


V G 5 -^ 0, vn = 0 = 


a 

Via = u. 


-i A -X 

It maj'' be noted that by setting u = 0, x^re make ^ = p for 
convenience. Under these conditionSy the kinetic equation for 
electrons can be written as, frem Eq. (4. IS) 2 


SZ + L.?4 = II ^ ve) ^ 

9>t ^ m Fq 3G t '‘m ^ 80 

0 © 


US. 

sp 


30 


VG 

© 


( 4 . 40 ) 

n p^ 
e-^ 


It is easily seen that if we multiply Eq. (4. 4jd hy Q - > 

G 

and integrate over dp, we obtain Eq. (3.11) for the conservation 
of energy. The integral on the right hand side of Eq. (4.4^1 
will drop to zero by virtue of Bq» (4.24). The x—c omp onent of 
the heat flux vector J is given by Eq^ (3.8) as j 
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2 -i 

■^ / p p F dp 

^ --CO 

e 


( 4.443 


Utilizing P from Eo,, (4.10), we obtain, 


J 


^ IT- 
6 e 9 6 


O rv,2ri 9 x 

2 m 6 
e 


CO ^ 

/ p p F^ hg (p) dp 


(4.4:^ 


By using the Maxwellian distribution for F'^, and substituting 

pi = _£ — ag Sections (4.5.1) and (4.5.2), we get i 

Vm© 


J 


X 


% ''V_j8 8_e 

(2iT)^72 (2 3 x 


where 


I 


© 


= f (P’)^ p' 



h^ (p*) dp' 


(4.4#) 


(4. 4-'!^) 


Using X from Eq. (4.8) and comparing Eq. (4.4-^ with the 
rate Eq. (3.12), we obtain the desired expression for thermal 
conductivity, K; 


9 "s k &V2 

2i% (Z+l) Unh m^ 


(4.4<b': 


This is found to be in agreement with the results of Spltzer 
and Harm [ 21,22] based on Fokker-Planck theory and also in agreement 
with the dominant term in the expression obtained by Sundaresan 
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and ¥u [20] "by a first order perturbation solution of BBGKf 
equations. 

The next chapter summarizes the present work reiterating 
the results and conclusions. 
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Chapter - 5 


SUMHARY^ 


C0NCLUUI01J3 and RECOMi 


IjpIjja r Julies 


5 » ^ pp e Corciu r,-:^-. 

In this work, transport coefficients of external 
fieo-d— free fully ionized "Dlasinas have h'.sn investigated. 

The transport coefficients studied are diffusion constant, 
viscosity and thermal conductivity. The study is based on 
a kinetic description of plasma. The results obtained are 
found to be in good agreement with those obtained by other 
authors based on different approaches, mainly F okker-Planck 
approach. 

In Chapter 1, the various models available for the 
analysis of plasma behaviour were review? ed, and a summary of 
background literature was presented. The phenomenological 
fluid approach and the kinetic approach were compared in 
some detail. The inadequacy of the Boltzmann kinetic equation 
for plasmas was pointed out. This is maijily due to long-range 
nature of the Coulomb force encountered in plasmas. 

The alternative kinetic approach, based on BBGEJ 
hierarchy of equations derived from lioville equation, is 
then presented in detail in Chapter 2. The problems of multiple 
collisions, and inter-particle correlations (statistical 
dependence), both arising mainly due to long range of the 
Coulomb force, were discussed. Finally Bogoliubcv’s approach 
to the solution of BBSKI equations was discussed. 
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In Chapter 3, appropriate moments of the kinetic equation 
are obtained leading to the definition of macroscopic variables, 
flux vectors and the conservation equations. The rate equations 
are introduced and transport coefficients defined. 

Chapter 4 presents our work. In this chapter, a simple 
method to obtain expressions for transport coefficients is 
developed. The method starts with the recognition of following 
facts; 

(1) The concept of transport coefficients is valid for 
n6;'r-- equilibrium situations, i.e., for systems in which deviation 
from equilibrium is small. Furthermore, for the concept to be 

of any practical use, the transport coefficients must not depend 
on the details of the deviation. Also due to the smallness of 
the deviation, only the first power of the magnitude of deviation 
need be considered. 

(2) Since the details of the deviation frcm equilibriiam 
are not important, the final expressions for the transport 
coefficients are essentially a result of equilibrium distribution 
functions. 

This is brought out in the present work. We start 

with small deviations from one-particle and two-particle distri- 
bution functions [Eqs. 4.10 and 4.14]. These when substituted 
into the first of BBGKI hierarchy of equations [Sq. 2.23], lead 
to the generalized kinetic equation [Eq. 4.19] which is no longer 
coupled with the remaining equations in the hierarchy. It is 
then possible to take the moments of the kinetic equation, once 
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again leaving out tlie details of the na-ture of deviation from 
equilibrium, represented by the functions h^ (p) in Eq. 4.10. 
Comparison with the rate equation then immediately leads to 
appropriate expressions for the transport coefficients. These 
coefficients as obtained in this work [Ens. 4,30 , 4.39, and 
4.45] can be summarized as followst 


1. Diffusion constant, D 


D oC 


Q 


5/2 




2. Viscosity, [j. 


|i oC 


\\ e 


5/2 


3. Thermal Conductivity, K 


K o( 


k e 


5/2 


V“e ® 


(5.1) 


(5.2) 


(5.3) 


These expressions are in agreement with the work of other 
authors based on P okker-Planck equation [21, 22, ’9]. -In 
case of thermal conductivity, Sundaresan and ¥u [20] have 
derived an expression foil owing Bogoliubov ' s approach to 
the solution of BBGKT equations. The expression obtained 
in this work agrees with the dominant term in their expression. 
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The constants of proportionality in Eqs. (5.1) to 
(5.5) can, in principle, be determined by solving the kinetic 
equation for hg(p) and then determining the integrals I , I 
and in Eqs. (4.30), (4. 35) and (4.45), respectively. This, 
however, is very lengthy and mathematically involved task, 
as is indicated by the work of Sundaresan and Wu [20], It 
should be noted that once the expressions in Eqs.(5.1) to (5.3) 
have been obtained, the constant of proportionality can be 
determined by a single experimental measurement in each case. 
Reference to Eqs. (4.30), (4.3^) and (4.45) indicates that 
the constants of proportionality will, in general, be functions 
of Z. This should be kept in mind if the same are to be obtained 
exp er im en tally . 

5.2 Suggestions for Future l/orh 

Instead of stiggesting minor variations of what has been 
reported here, we will end this work by pointing out areas in 
which further progress will be of maximum value % 

1. For plasmas with large gradients in the macroscopic 
quantities, for example, the plasma on the surface of a laser 
irradiated fusion-pellet, the hydrodynamic approach based on 
the first five moments of the distribution function, namely 
n, u and T, is not valid. Hence a method based on either 
additi cnal moments or a direct solution of the initial value 
kinetic problan should be attanpted. It should be noted that 
the existing hydrodynamic codes for laser— pellet— fusion circumvent 
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this difficulty by empirical concepts like heat- inhibit ion 
factors etc. This appears to be unsatisfactory. 

2. The classical description of plasma does not seem 
to be satisfactory due to divergence x)-S'0’blem at small inter- 
particle distances. Quantum effects should be included. 

3. Methods should be extended to partially ionized 
plasmas, with man^ species. This will necessitate the modelling 
of ionization-recombination behavior ana charged-neutral as well 
as neutral-neutral collisions, which can be a very complicated 
task to perform satisfactorily. 

4. External fields, either in the form of steady 
electric and magnetic fields or electromagnetic radiation, or 
something in-between, should be considered. 

5. Relativistic effects in the plasma, necessitating the 
consideration of interparticle magnetic interaction, the plasma 
radiation, and itvS interaction with plasma particles, should be _ 
included. 

6. It is impossible to keep a plasma free from collective 
phenom ena( linear and nonlinear waves). In fact the collective 
mechanisms can be dominant modes for transfer of mass, mcmentum 
and energy, and can manifest in the form instabilities. The 
list of collective modes that can be excited in a plasma is very 
large. Methods to specify the collective modes that are likely 

to be excited under given conditions and their effect on transport 
phenomena in a plasma should be included in ary realistic analysis- 
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7. The pi’Oblems related to turbulence and fluctuations 
in plasmas have been attacked only recently and there is 
sufficient scope for further work. 

8. There are problans of more basic nature, which can 
be attacked. These relate to the construction of an irrever- 
sible statistical or kinetic theory based on a reversible 
dynamic theory, insufficiency of the statistical element in the 
distribution fimctions, etc. These problems need to be investi- 
gated further. 

9. Pinaliy, any theory for denso plasmas is likely 
to be very complica.ted. It will, therefore, be necessary to 
supplement on substantiate theoretical effort by experimental 
work aimed at measuring plasma properties under various 


conditions. 
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Appendix I 

Properties of Equilibrium Distribution Function 
1 . One~P article Equilibrium Distribution Enaction 

■where 

^ = p - mu, .r = 1 . 

The following properties of E® are easily verified; 


00 


/ E^ dp = 1, 

oo 

(I"l) 

00 

/ dp = 0, 

- CO 

(1-2) 

00 Q 

= 1®. 

(1-5) 

1 

II 

O 

(1-4) 

CO 

/ dp = 0, etc, 

-oo 

(1-5) 


In tk ; above integrals, it may be noted tliat 

= dp^ dpy dp^ = dp = 

The properties (I-l) to (1-5) can be verified, either in the 
Cartesian coordinates or by switching over* io spherical polar 


coordinates. 



By definition, we must have 
/ ^P2 ‘^^2 VB(1) 

It is easily seen that the equilihrixun two-particle 
distribution function given by Eq. (2*59) satisfies th 
above identify if E > > 1, which will always be the 


case. 
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Appendix II 


lay Ion Series Expansion of C(r, 0 , nJ 


let 



a, then, 


G(r, 6-, n) = 1 


0r 


exp(- ") , 

l-r^ 


where 

e 9 »/2 

r = ( -^ ) (see Eq. 2.31) • 

ne^ 


(II. 1) 


By direct partial differentiation, we obtain, 


^ _ ^^1^2 r 

59 - q2^ L -L 



] exp (- 



and 


c>n 


aZ^Zj 2 
20rj^n' 


exp ( 



o 


(II. 2) 


(II .5) 


Now, we write, 


G(r5 9, n)i!i G(r , 





(e - 9;^) 

(n - n^) , 

(II .4) 


e-e 


®2“®1 


1 


*7 G»r 


1 dQ 
2^ 


X, 


where 


(II .5) 
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and 


n-n^ = 


2 



n 


A 

■r 


1 ^ X 

2 S3c 


(n.6) 


•where x denotes the x-component of r. 

Substituting Eqs. (II. 2), (II. 5), (II.5) and (II 6) into 
Eq. (II. 4) > we get the desired expansion for Gr(r, 6, n) as 
follows; 


G(r, ©•, n) = 


aZiZp 

1 " exp( - — ) 

1 ^D1 


, n r . . , r ^ ^ © I 

"lif 


X 

r 


-h 


4© 


^D1 1^ 


X 


■D1 


(II. 7) 


where, 


•El 


( ) 

n^e Z 


1/2 


Thus G(r, §, n) has been expressed in terms of -the 
macroscopic variables at point 1 and their derivatives at 
the same point. 
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Appendix III 


Evaluation of the Collision Integral in Eg, (At,17 ) 


In this appendix we v/ill show that Eq. (4.17), nanely 


I = n 


“ _i ^ 0 — A 

/ dq_ G(r, e, n) / dpp F(2), 

9q^ - CO ^ 

(III.l) 


reduces to 


T - 9 F r a. Zii 1 

39 ^^'LQ+n-l* 

9p 


By definition, the integral J dpp F(2) is unity. To evaluat- 

^00 ^ 

the other integral, we note that the interparticle potential, 
as given by Bq. (4.16) is 


^2 ® 


so that 


1 ir 


(III. 2) 


wJI 

91 ■ '92 


+ (<liy-42y)3 + 


(<1ix-92e)^+ ■*■ (iiwW > 
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Fow, from Appendix II we have 

8.Z1-1 Z , 


G(r, 0, H) = [1 - gxp(~ 


•D1 



+ 

2 

1 1 

H 

1 

r 


2 

29^ 

2.x 

^ D1 


+ 

aZ,Z^ 

exp(. 

r ■ 

^Dl 

can be 

writt ePx 

as 



G(r, 

II 

1 

iCD 

G1 + 

G2 + 

G3, 


3X 


T_) ^ 
-Dl' 


X 


, X 


1 


1 "’Dl 


where G1 is the first bracketed term, G2 and G3 the second 
and third terms, respectively. 

Then ef ore, 


.-rp/T-v Z-iZo© “ V 

I = _ n —I—I [ / dq^ G1 + / -3 G2 

3 Pi o -® r^ r 


+ / 4o ^ e3]. 


^2 3 

r 


(III. 4) 


The first integral of Eq. (III. 4) is, 


CO ^ X* 


(III. 5) 



z 


6: 


S inc e r = 


dr = 


X a 

" ^2 


- dq. 


r->+“ as q 2 
so the expression (III, 5) takes the form 


yl 



Ftu Hl-j 


“ S-Z-i Zo r> , ^ 

/ -^ [1 - exp(- dr. 


[1 


•D1 


(III. 6) 


Transforming to spherical coordinates as shown in Pig. I'^'I 

¥0 have? 


s 

r 


.[ i r sin0 cos i|^ + 3 r sin0 sin 'i' + k cos0], 


and dr = r^ sin0 d0 d^* dr 


In view of the above substitutions, the integral becomes, 


% 


2r 


aZf-i 2 j r\ Y* 

/ ^ [1 --557 


)]dr / sin^0d0 / cos d^- 


% 2% 
2r 


- f J [1 


^^1^2 (_«£-) ]dr / sin^0 d0 / sin ^ dtJ) 


- ^ exp, _ 

^D1 


IX 


2 jx 


j. 

/ k [1 


aZ' 


1^2 /„.X-)]dr J sin0cos(/d(2f dll'. 

^Bl - 

o o 



> I « ■ ai I 


I . m 
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By direct integration it can be shown that, 

2 n. 2k Tt 

I cos4'dii^ =0, / sin i!) dll' = 0 and / s±n<^ cos0 d<^ = 0 . 

0 O o 


Therefore, the first integral in Eq. (III. 4), that is 


J ^ G1 dig, 

- CO ^ 


drops to zero. 

The second integral is 


jk 

r 

"5 


G2 dr 




2 © 


2 3 x 


c» ^ 

/ ^ 

l~oo p- 


[1 


2r. 


D1 


'] exp(- 


■D1 


X 

r” 


dr 


[/ i ( i . 

2©2 ^^^1 ^ 


2r. 


D 1 


•) exp(- — )dr' 


^ D 1 


It 2tt 

/ sin^(2f d(? I cos d>jj 


0 


o 


7U 2tc 

,3. 


+ r i*(i_ — ) exp(- - — )dr / sin^(2^dJ5 / simile os W 

0 


o 


2 % 


2k 


A 

+ / k 
o 


(1- p^) exp(- f--)dr / sin^gcos^d^i’cos^d^'l 
D1 D1 

0 o 


(III. 7) 
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Evaluating the integrals it can he shown that 

2n 

f sin 4* c os d = 0 and / c os 'I' d ’f' 


0 


o 


n 2% 

and to find / sin^^ d0 / cos ’-1^ d ^ , -we' make use of 
o 0 

the fact that 


sin^0 = I sin0 - j sin50, 


so that. 


n 


/ sin^0 d0 


4 

3 ’ 


and since 


2jz 


2 % 


/ cos^ rp dilj = / (1 + cos2 4) d^ = 


7C 


¥e have. 


n 2 k 


/ sin^0 d0 / cos^t}) d4 
o 0 


■ ~4k 

3 


With these results Eq. (III. 7) takes the form 

fr> 

^^i£ 2 _© 1 A I /-i _r N r n 


( ^ G2 dr = - ^ ^ 

3 2 X j 3 


.[V- 


2r. 


D1 


■) exp(- - — ) dr 


■D1 
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ZaZ^ZgH ag 
2 ^3: 

3 ©1 


± 


if- 

a 


-r/r. 


D1 


dr 


I f ^ 
- 2j 
<> 


e ^ dr], 

^D1 


Let us put 


z; s, 


-D1 


then the above equation takes the form 


2aZn Z nit 


lisi i® I 

^ ©2 b x|^ ^ 


9 

j 


e-= ds - -f- 


I ^ 


s e 


ds] 


2az^a_^ 


3 e 


1 


C?x I 


3 ( 2)1 


since T (1) = p (2). 




3 © 


0x L ^ • 


1 


Hence we have 


f 


T» 

— ^ Gn = 


- t 


•Pl 

3 


aZ^Z2it 


a© 

c>X 


(III. 8) 


Finally we need to evaluate 




1 ^P1 


dr 
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TC 


2n 


in 


I ^ J e 2 :p(-- ^ — )dr j" sin^{2fd0 j" cos^-^ d f 


% 


O 

2k 


+ 


5 - I r exp(- — ) dr 1 sin^^ d0 j sin 't cos 4 d*f 
i ^Dl 


> 


n 

.. 2 r. 


2k 


+ k \ r exp( ) dr \ sin^0 cos0 d0 

4 ^D1 J 


J c os 4 d 4 ] • 


The last two integrals vanish and we have 


^^1^2^ ±n 1 ^ 




1 [ J r exp(- ^ — ) dr] 


D1 


= _ 5: ^ZiZgTtrpi ^ I 

3 ^ h 


So 
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G_ dr 
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A 

- X 


a Zj^Zgnrpi | 


(III. 9) 


Putting Bqs, (III, 8) and (III,9) into Bq. (III. 4) we finally 
get 
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nZ-, Z^e 


l'^2^ ^B(l) r - f D1 ^^1^2^ <1© I 
^^o '‘- 3 a2 ) 


Pi 


r^^T aZj^Z^-Ji 

91 


^D1 s-Zj^Z^-n: 
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2 2 2 
na 

3-0 




'\i e 
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4 * 


^ n - 
n - 


(III. 10) 


where we have dropped the subscript 1, and used the 
not ation. 


a 



(as defined in Appendix II). 


The ri^t had number of Eq. (II 1. 10) is the evaluated 
expression for the collision integral appearing in Eq. (4.17). 
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